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ABSTRACT 


A  family  of  self-similar  fields  provides  the  two  parameters  required  to  characterize 
the  full  range  of  high-  and  low-triaxiality  crack  tip  states.  The  two  parameters,  J  emd  Q, 
have  distinct  roles;  J  sets  the  size  scale  of  the  process  zone  over  which  large  stresses  and 
strains  develop,  while  Q  scales  the  near-tip  stress  distribution  relative  to  a  high  triaxiality 
reference  stress  state.  An  immediate  consequence  of  the  theory  is  this:  it  is  the  toughness 
values  over  a  range  of  crack  tip  constraint  that  fully  chEU-acterize  the  materisd’s  fracture 
resistance.  It  is  shown  that  Q  provides  a  common  scEile  for  interpreting  cleavage  fracture 
and  ductile  tearing  data  thus  allowing  both  failure  modes  to  be  incorporated  in  a  single 
toughness  locus. 

The  evolution  of  Q,  as  plasticity  progresses  from  small  scale  yielding  to  fully  yielded 
conditions,  has  been  quantified  for  several  crack  geometries  and  for  a  wide  range  of  material 
strEun  hardening  properties.  An  indicator  of  the  robustness  of  the  J-Q  fields  is  introduced; 
Q  as  a  field  parameter  and  as  a  pointwise  measure  of  stress  level  is  discussed. 
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1.  INTRODUCTION 


A  two- parameter  fracture  theory  can  be  motivated  by  considering  the  progression  of 
plastic  states  zis  loading  on  a  cracked  body  is  increased.  At  low  loads,  the  near-tip  stresses 
and  deformations  evolve  according  to  a  self-similar  field,  scaled  by  Rice’s  J-integral  [1]. 
This  field,  characterized  by  a  high  level  of  stress  triaxieJity,  also  describes  the  evolution  of 
the  near-tip  stresses  zind  deformations  in  certain  crack  geometries  as  plastic  flow  progresses 
from  well- contained  yielding  to  Izu'ge  scale  yielding.  While  this  high  triaxiality  field  is 
only  one  of  many  possible  states  that  can  exist  under  fully  yielded  conditions,  it  is  the 
only  field  that  has  received  careful  study  until  recently.  When  the  high  triaxiality  field 
(Hutchinson  [2],  Rice  and  Rosengren  [3],  Rice  and  Johnson  [4],  McMeeking  [5])  prevails 
over  distances  comparable  to  several  crack  tip  openings,  J  alone  sets  the  near-tip  stress 
level  and  the  size  scale  of  the  zone  of  high  stresses  and  deformations.  Considerable  efforts 
have  been  directed  to  establishing,  for  different  crack  geometries,  the  remote  deformation 
levels  which  ensure  that  the  near-tip  behavior  is  uniquely  measured  by  J  (McMeeking  and 
Parks  [6],  Shih  and  German  [7]).  The  end  result  is  a  framework,  based  on  J  and  the  high 
triaxiality  crack  tip  field,  for  correlating  crack  growth  over  a  range  of  pleme  strain  yielding 
conditions  (see  review  articles  by  Hutchinson  [8j,  Parks  [9])  and  for  relating  critical  values 
of  the  macroscopic  parameter  Jjc  to  fracture  mechanisms  operative  on  the  microscale  (see 
review  article  by  Ritchie  and  Thompson  [10]). 

Arguments  that  a  single  parameter  might  not  suffice  to  cheuracterize  the  near-tip  states 
of  fully  yielded  crack  geometries  have  been  raised  by  McClintock  [llj.  He  noted  that 
non-hardening  plane  strain  crack  tip  fields  of  fully  yielded  bodies  are  not  unique  but 
exhibit  levels  of  stress  trizodality  that  depend  on  crack  geometry.  Though  high  stress 
triaxiality  is  maintained  in  geometries  involving  predominantly  bending  on  the  uncracked 
ligament,  the  level  of  crack  tip  stress  triaxiality  in  geometries  dominated  by  tension  loads 
generally  decreases  as  yielding  progresses  into  the  fully  plastic  state  (see  Refs  6  and  7). 
Indeed,  experimentally  measured  7— resistamce  curves  for  center-cracked  panels  exhibit 
significantly  higher  slopes  than  those  for  compact  specimens  (Begley  and  Landes  [12]). 

Undoubtedly,  there  is  merit  to  the  latter  point  of  view  and  it  is  the  purpose  of  this 
paper  to  show  that  this  viewpoint  can  be  properly  reconciled  by  a  two-parameter  J  —  Q 
theory  proposed  by  O’Dowd  and  Shih  [13,14]  and  Shih  ct  al.  [15,16].  They  showed, 
through  full-field  finite  element  calculations,  that  the  J  —  Q  fields  dominate  over  physically 
significant  size  scales,  i.e.  they  represent  the  environment  in  which  the  ductile  and  brittle 
failure  mechanisms  are  operative.  An  approach  based  on  higher-order  asymptotics  has 
been  suggested  by  Li  and  Wang  [17]  and  pursued  by  Sharma  and  Aravas  [18]  and  Chao 
tt  al.  [19].  Extending  the  analysis  in  Refs  17  and  18,  Xia  et  al.  [20]  have  obtained  up  to 


1 


five  terms  in  the  asymptotic  series  and  showed  that  the  collective  behavior  of  the  series  is 
consistent  with  the  J  —  Q  field. 

An  alternative  two-peirEimeter  approach  based  on  J  and  the  elastic  T— stress  has  been 
advocated  by  Betegon  and  Hancock  [21],  Al-Ani  and  Hancock  [22],  Du  and  Hancock  [23], 
Parks  [24],  Hancock  et  al.  [25]  and  Wang  [26,27].  Under  circumstances  where  it  is  applica¬ 
ble,  the  J  —  T  theory  can  be  shown  to  be  equivalent  to  the  J  —  Q  theory.  This  is  discussed 
in  the  section  on  small  scale  yielding.  The  toughness  scaling  approach  of  Dodds  et  al.  [29] 
can  also  be  shown  to  be  consistent  with  the  J  —  Q  theory  (see  Kirk  et  al.  [30]). 

2.  J-Q  THEORY 

Consider  a  cracked  body  of  characteristic  dimension  L  loaded  remotely  by  a  stress 
denoted  by  (7°°.  The  scale  of  crack  tip  deformation  is  measured  by  J/cxq  where  <7o  is 
the  material’s  tensile  yield  stress.  It  can  be  shown  from  dimensional  grounds  that  when 
L  ^  J/(Tq  all  neax-tip  fields  axe  members  of  a  single  family  of  crack  tip  fields.  Each  member 
field  is  characterized  by  its  level  of  deformation  as  measured  by  J/<tq  and  its  level  of  crack 
tip  stress  triaxiality,  as  measured  by  Q,  which  also  identifies  that  field  as  a  particular 
member  of  the  family.  For  example,  the  self-similar  solution  of  Rice  and  Johnson  [4]  and 
McMeeking  [5]  or  the  HRR  field  (Refs  2  and  3)  can  be  taken  as  the  Q  =  0  member  field. 
In  short,  the  Q— family  of  fields  provides  the  proper  characterizing  parameters  for  the  full 
range  of  neax-tip  fracture  states. 

The  weak  coupling  between  deformation  and  stress  triaxiality  in  a  plastically  deform¬ 
ing  material  provides  another  argument  in  favor  of  a  two  peirameter  description  of  near 
tip  states.  Since  plastic  flow  is  incompressible,  the  superposition  of  a  purely  hydrostatic 
stress  state  induces  only  an  elastic  volume  change.  Consider  a  plastically  deforming  ma¬ 
terial  element  in  the  forward  sector  of  a  crack  as  depicted  in  Fig.  1.  We  can  superpose  a 
hydrostatic  stress  Qctq  with  little  or  no  effect  on  the  deformation  state.  It  follows  that  near 
tip  deformation  and  stress  triaxiality  cannot  be  scaled  by  a  single  parameter  such  as  J. 
A  second  parameter  is  required  to  quantify  the  level  of  crack  tip  stress  triaxiality.  Clearly 
this  argument  does  not  apply  to  the  back  sector  since  traction  free  conditions  must  be 
satisfied  on  the  crack  faces.  However  this  is  of  no  physical  consequence  since  the  fracture 
processes  occur  in  the  forward  sector,  which  is  therefore  the  region  of  interest. 

A  size  scale  must  enter  into  the  fracture  description.  We  focus  on  fields  ahead  of 
the  crack  that  aie  relevant  on  the  scale  of  the  crack  opening  displacement  St,  or  J/cto, 
representing  the  environment  in  which  the  failure  mechanisms  are  operative. 
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Back  Sector 


Forward  Sector 


f  <3cro 


Figure  1.  Schematic  illustrating  the  necessity  for  a  hydrostatic  stress  parameter  zmd  a  de¬ 
formation  parameter  to  characterize  the  full  range  of  near  tip  states  in  the  forward 
sector. 

2.1  Q-Family  of  Fields  —  MBL  formulation 

It  proves  convenient  to  construct  the  Q— family  of  fields  using  a  modified  boundary 
layer  (MBL)  formulation  in  which  the  remote  tractions  are  given  by  the  first  two  terms  of 
the  small-displacement-gradient  linear  elastic  solution  (Williams  [31]), 

I,,- = +  TMiy.  (2.1) 

VZTrr 

Here  Sij  is  the  Kronecker  delta;  r  and  9  are  polar  coordinates  centered  at  the  crack  tip 
with  6  =  0  corresponding  to  a  line  ahead  of  the  crack  as  shown  in  the  insert  to  Fig.  2. 

Fielas  of  different  crack  tip  stress  triaxialities  can  be  induced  by  applying  different 
levels  of  r/cTo-  From  dimensional  considerations,  these  fields  cam  be  organized  into  a  family 
of  crack  tip  fields  of  the  form: 

where  J  is  Rice’s  J— integral  [1].  That  is,  the  load  parameter  T/cro  provides  a  conve¬ 
nient  means  to  investigate  and  parameterize  specimen  geometry  effects  on  near- tip  stress 
triaxiality  under  conditions  of  well-contained  yielding.  Indeed,  such  studies  have  been 
carried  out  by  Betegon  amd  Hancock  [21],  Bilby  ei  al.  [32]  and  Haxlin  and  Willis  [33]. 
Nevertheless,  the  result  in  Equation  2.2  cannot  have  general  applicability  since  the  elaistic 


3 


solution  (Eq  2.1)  upon  which  the  T— stress  is  defined,  is  an  asymptotic  condition  which  is 
increasingly  violated  as  plastic  flow  progresses  beyond  well-contained  yielding. 

Recognizing  the  above  limitation  O’Dowd  and  Shih  [13,14],  henceforth  referred  to  as 
OS,  identified  members  of  the  plane  strain  family  of  fields  by  the  parameter  Q  which  arises 
naturally  in  the  plasticity  analysis.  OS  write: 


The  additional  dependence  of  fij,  gij  and  hi  on  dimensionless  combinations  of  materizJ 
parameters  is  understood.  The  form  in  Equation  2.3  constitutes  a  one-pzirameter  family 
of  self-similar  solutions,  or  in  short  a  Q-family  of  solutions.  The  annular  zone  over  which 
Equation  2.3  accurately  quantifies  the  actual  field  is  called  the  J  —  Q  annulus.  Represen¬ 
tative  distributions  of  the  Q  family  of  fields  are  presented  in  Fig.  4  of  Ref  1  and  Fig.  1  of 
Ref  2. 

2.2  Difiference  Field  and  Near-Tip  Stress  Triaxiality 

Using  the  modified  boundary  layer  formulation,  and  considering  a  piecewise  power 
law  hardening  material,  OS  generated  the  full  range  of  small  scale  yielding  plane  strain 
solutions,  designated  by  ((Tij)ssY-  OS  considered  the  difference  field  defined  by 

icrij)dia  =  (<7’,y)ssy  -  (<^«i)HRR,  (2-4) 

where  (c7,;)hrr  is  the  HRR  field.  They  systematically  investigated  the  difference  field 
within  the  forward  sector,  |^|  <  7r/2,  of  the  annulus  J/ctq  <  r  <  5J/<7o,  since  this  zone 
encompasses  the  microstructurally  significant  length  scales  for  both  brittle  and  ductile 
fracture  (see  Ref  10).  Remarkably,  the  difference  field  in  the  forward  sector  displayed 
minimal  dependence  on  r.  Noting  this  behavior,  OS  expressed  the  difference  field  within 
the  forward  sector  in  the  form 


(2-5) 

where  the  angular  functions  Cij  are  normalized  by  requiring  d^e${6  =  0)  to  equal  unity. 
Moreover,  the  angular  functions  within  the  forward  sector  exhibit  these  features:  ~ 

ss  constant  and  \or0\  <  \^e9\  (see  Figs.  3,  4  and  5  in  Ref  13). 

Thus  the  difference  field  within  the  sector,  1^|  <  7r/2  and  J/ctq  <  r  <  5J/ao,  corre¬ 
spond  effectively  to  a  spatially  uniform  hydrostatic  stress  state  of  adjustable  magnitude, 
i.e.  ((7,j)djff  =  QcoSij.  Therefore,  Q  defined  by 

^  ^  at  ^  =  0,  r  =  27/<to  (2.6) 

(To 
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is  a  natural  measure  of  near-tip  stress  triaxiality,  or  crack  tip  constrmnt,  relative  to  a  high 
triaxiality  reference  stress  state.  In  words,  Q  is  the  difference  between  the  actual  hoop 
stress  and  the  corresponding  HRR  stress  component,  the  difference  being  normalized  by 
Cq.  For  definiteness  we  have  evaluated  Q  at  r  =  2J/cro,  how'ever  w'e  point  out  that  Q  is 
effectively  independent  of  distance.  The  distance  chosen  for  the  definition  of  Q  lies  just 
outside  the  finite  strain  blunting  zone  so  that  Q  from  a  small  or  finite  strain  analysis  should 
be  nearly  the  same. 

OS  also  considered  the  difference  field  whereby  the  standard  plan<  strain  small  scale 
yielding  solution  (o’,j)ssy;T=0j  which  is  driven  by  K  alone,  is  the  reference  solution,  i.e. 

(«^o)diff  =  i<^ij)sSY  -  (<7’,j)sSY;T=0-  (2-7) 

In  this  case  the  difference  field  in  the  forward  sector  matches  a  spatially  uniform  hydrostatic 
stress  state  even  more  closely.  Thus  an  alternative  definition  of  Q  is 

at  «  =  0,  r  =  2J/<ta.  (2.8) 

(To 

A  definition  of  Q  consistent  with  its  interpretation  as  a  triaxiality  parameter  is, 

^<Tm]_SSY;T=0  ^  ^  ^  ^  (2.9) 

(To 

where  <Tm  is  the  hydrostatic  stress.  We  have  calculated  Q  based  on  the  hoop  stress  (Equa¬ 
tion  2.8)  and  the  mean  stress  (Equation  2.9)  for  the  full  range  of  T-stresses  and  finite 
width  geometries.  We  have  found  that  the  difference  between  Q  and  Qm  is  always  less 
than  0.1.  While  the  values  of  Q  presented  in  this  paper  are  ceJculated  from  the  hoop  stress 
via  Equation  2.8  it  is  clear  from  the  above  that  these  Q  values  can  also  be  used  to  calculate 
the  corresponding  hydrostatic  stress  levels. 


2.3  Difference  Field  and  Higher-Order  Terms  of  the  Asymptotic  Series 

The  connection  between  the  difference  field  and  higher-order  terms  of  th*  as\  .nptotic 
series  can  be  understood  in  the  context  of  the  MBL  formulation.  Here  the  stress  lield  obey 
the  functional  form 

(T.y  =  (Tofij  ( ’  (2-10) 

which  also  should  apply  to  finite- width  crack  geometries  as  long  as  the  characteristic  crack 
dimension  L  is  sufficiently  large  compared  to  J/cro.  Now,  if  one  assumes  a  product  de¬ 
pendence  on  the  first  argument  in  Equation  2.10  and  works  within  deformation  plasticity 
theory  and  an  elastic  power-law  hardening  material,  then  one  obtains  a  series  in  r/  ( J/ao ): 
/  j  \i/(n-(-i) 

(Tti  =  (To  I  - z —  I  d’i,(^;  n)  -I-  second  order  term  -I-  higher  order  terms  (2.11) 

\Qeo<roInrJ  " - ' 

Difference  Field 
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where  cq  is  a  reference  strain,  a  a  material  constant  (for  a  piecew'ise  power-law  material 
Q  =  1)  cind  In  is  an  integration  constant.  By  definition,  the  asymptotic  series  beyond  the 
first  term  is  equivalent  to  the  difference  field  since  (see  previous  section) 

i7,j  =  (cT,j)HRR  +  Difference  Field.  (2.12) 

The  HRR  field  and  the  second-order  term  provides  only  a  two-term  approximation  to  the 
solution  for  the  MBL  problem  and  this  point  appears  not  always  to  be  understood. 

The  higher-order  asymptotic  cinaJysis  of  Refs  17  and  18  has  been  extended  by  Xia, 
Wang  and  Shih  [20].  They  have  obtained  a  five  term  exp2insion  for  the  series  in  Equa¬ 
tion  2.11  for  n  =  3  and  four  term  expansions  for  n  =  5,7  and  10.  Furthermore,  they 
have  successfully  matched  the  four  term  series  to  the  radial  and  angular  variations  of  the 
difference  field  given  in  Fig.  3  and  Fig.  5  in  Ref  13  for  an  n  =  10  material.  Indeed  in 
the  forward  sector  |5|  <  7r/2,  the  collective  behavior  of  the  the  second,  third  £ind  fourth 
order  terms  is  effectively  equi%’alent  to  a  spatially  uniform  hydrostatic  stress  state.  This 
obserN'ation  together  with  the  discussion  of  the  previous  section  supports  the  simpler  form 
for  the  near-tip  fields: 


^0  =  (<^u)hrr  +  QtroSij,  1^1  <  7r/2.  (2.13) 

Furthermore,  it  may  be  noted  that  an  admissible  range  of  stress  states  for  an  elastic- 
perfectly  plastic  material  can  be  written  in  the  form 

(^ij  =  (c^.i)PrandU  +  Q<TQ6ij  ,  |5|  <  7r/4,  (2.14) 

DifTercDce  Field 

where  (ffijOprandtl  designates  the  Prandtl  slip-line  solution  and  again  the  difference  field 
corresponds  simply  to  a  uniform  hydrostatic  stress  state  scaled  by  Q  (see  Refs  23  and  14). 

2.4  Variation  of  Q  with  Distance 

Because  Q  scales  the  difference  field  relative  to  a  reference  stress  state,  it  provides 
a  sensitive  measure  of  the  evolution  of  near-tip  stress  triajdaJity  in  finite  width  cracked 
bodies.  It  also  can  be  used  to  detect  changes  in  the  stress  triaxiality  that  deviates  from 
the  pattern  that  develops  under  MBL  loadings.  For  this  purpose,  we  consider  Q(f)  defined 
by 

Q(f )  ^  <’••(^'1  -  l'^«»(^)lssY.T..  ^  ^^0  (2.15) 

(To 

where  f  =  r/{J/(7o)-  Note  that  (cr«« )ssY;T=o  is  chosen  as  the  reference  field. 
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The  mean  gradient  of  Q  over  1  <  r  <  5, 


Q'  = 


Q{f  =  5)-Qif  =  l) 
4 


(2.16) 


can  be  used  to  monitor  changes  in  the  spatial  pattern  of  the  stress  triaxiality  ahead  of  the 
crack  that  do  not  conform  to  a  spatially  uniform  hydrostatic  stress  field.  In  other  words, 
Q'  provides  a  measure  of  the  robustness  of  the  J  —  Q  fields  for  the  application  of  interest. 
If  IQ'l  <  0.1  then  the  difference  field  is  effectively  constant  over  1  <  f  <  5.  If  |Q'|  is  much 
larger  them  0.1  then  the  variation  of  the  difference  field  over  the  interv’zJ  1  <  f  <  5  cein 
be  comparable  to  <tq.  In  the  latter  case,  Q  provides  only  a  pointwise  measure  of  the  stress 
level — at  the  distance  2J/<tq — ahead  of  the  crack  tip. 


2.5  Reference  Field  Distributions 

Table  1  provides  the  reference  field  distributions  for  a  broad  range  of  n  values.  For 
completeness  we  have  included  the  hoop  stress  distributions  according  to  the  HRR  singu¬ 
larity  and  the  small  scale  yielding  solutions  for  small  straiin  and  finite  strain  (piecewise- 
power-law  hardening  material;  EI<tq  =  500,  i/  =  0.3).  Figure  2  presents  the  hoop  stress 
and  mean  stress  reference  fields  established  by  the  MBL  formtilation  with  T  =  0.  The 
originsJ  studies  of  OS  were  based  on  a  finite  strain  formulation  to  ensure  a  full  description 
of  the  near-tip  states.  Our  subsequent  studies  have  shown  that  small  and  finite  streun 
analyses  provide  essentially  identical  resvdts  over  the  region  of  interest  1  <  f  <  5.  This 
can  also  be  seen  by  comparing  the  finite  strain  and  small  strain  distributions  in  Fig.  2. 

Two  reference  fields  have  been  proposed, ,  (<rjs)HRR  and  (<7's^)ssY;T=o-  Our  numericeJ 
investigations  of  different  crack  geometries  show  that,  when  the  small  scale  yielding  solution 
is  chosen  as  the  reference  state,  the  difference  fields  correspond  more  closely  to  a  imiform 
hydrostatic  stress  state  over  a  greater  rwge  of  plastic  deformation.  However,  the  choice 
of  reference  distribution  used  in  the  definition  of  Q  remains  a  matter  of  convenience.  We 
emphasize  that  once  a  choice  is  made,  it  must  be  applied  consistently  throughout  the 
analysis.  Wc  recommend  that  Equation  2.8  be  used  as  the  standard  definition  for  Q  with 
the  small  strain  solution  as  the  reference  field.  Having  a  standard  definition  facilitates  the 
comparison  of  solutions  obtained  by  different  investigators  and  the  tabulation  of  a  handbook 
of  Q  solutions. 

While  we  have  limited  our  discussion  to  a  piecewise  power-law  hardening  material, 
the  J  —  Q  theory  is  independent  of  the  form  of  the  material’s  constitutive  relation.  Thus 
(<7^^ )ssY;T=o  cm  be  evaluated  for  an  actuzd  stress-strain  relation.  Of  course,  for  consistency 
the  malyses  in  the  fracture  application  should  also  use  the  same  stress-strmn  relation. 
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Table  1.  Reference  stress  distributions,  ceefaQ,  from  HRR  field  and  small  and  finite  strain 
boundary  layer  solutions. 
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2.6  Engineering  Applications  of  the  J  —  Q  Theory 

For  engineering  applications  two  forms  of  the  near  tip  pleistic  states  axe  proposed: 

=  (<^o)hrr  +  (2-17) 

and 

<7ij  =  (<^fi)ssY:T=o  +  (218) 

where  Q  in  Equations  2.17  and  2.18  are  defined  by  Equations  2.6  and  2.8  respectively. 

The  v’alues  of  the  hoop  stress  of  the  HRR  field  for  1  <  r  <  5  is  given  in  Table  1.  The 
other  stress  components  can  be  found  in  Symington  ct  al.  [39].  The  hoop  and  mean  stress 
distributions  of  the  small  scale  yielding  field  with  T  =  0  are  given  in  Table  1  and  Fig.  2. 
More  details  are  fotmd  in  Refs  13  and  14. 

The  physical  interpretation  of  Equations  2.17  and  2.18  is  this:  negative  (positive)  Q 
values  mean  that  the  hydrostatic  stress  ahead  of  the  crack  is  reduced  (increased)  by  Qgq 
from  the  J— dominant  stress  state,  or  the  standard  small  scale  yielding  stress  state.  This 
interpretation  is  precise  when  \Q'\  1. 

As  stated  previously  we  recommend  the  use  of  Equation  2.18  in  the  J  —  Q  fracture 
methodology.  However,  the  explicit  representation  in  Equation  2.17  cam  facilitate  analyses 
leading  to  predictions  of  constraint  effects  on  toughness  as  outlined  in  the  section  on  a 
cleavage  toughness  locus. 

3.  SMALL  SCALE  YIELDING 
3.1  Q-T  Relation 

Within  the  modified  boundary  layer  formulation,  a  strict  one-to-one  correspondence 
exists  between  Q  and  T  (see  Refs  13  and  14).  For  a  piecewise  power-law  hardening  material, 
the  relationship  takes  the  form 

(3  =  F(r/ao;n),  (3.1) 

with  an  additional  weak  dependence  on  E/o-q  and  I'j  where  E  is  Yoimg’s  modulus  and  v 
is  Poisson’s  ratio.  Curves  of  Q  vs.  T/<xq  for  n  =  3,  5,  10,  20  and  oo  axe  displayed  in  Fig. 
3.  The  Q  values  in  Fig.  2  and  Table  2  are  based  upon  the  definition  in  Equation  2.8, 
and  were  determined  by  small  strain  analyses,  using  E/cq  =  500  and  v  =  0.3;  essentially 
identical  results  were  obtained  from  finite  strain  analyses.  It  can  be  seen  that  Q  increases 
monotonically  with  Tja^.  Also  note  that  crack  tip  stress  triaxiality  can  be  significantly 
lower  than  the  reference  state  (the  Q  =  0  state)  but  cannot  be  elevated  much  above  it. 
The  values  Q  and  Q'  are  given  in  Table  2  keeping  only  two  places  beyond  the  decimal 
point.  Note  that  the  largest  value  of  Q'  is  less  than  0.04.  Thus  Q  is  effectively  constant 
over  the  distance  1  <  f  <  5  for  all  MBL  loadings. 
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TIco 

-1.0 

-0.75 

-0.5 

-0.25 

0.0 

0.5 

1.0 

n  =  3,Q 

-0.37 

-0.16 

0.0 

0.27 

0.41 

Q' 

O'* 

0.0 

0.0 

-0.01 

-0.04 

n  =  5,Q 

-0.74 

-0.46 

-0.20 

0.0 

0.27 

0.36 

Q' 

0.0 

0.0 

0.0 

0.0 

0.0 

-0.02 

n  =  10,  Q 

-1.23 

-0.92 

-0.54 

-0.23 

0.0 

0.21 

0.24 

Q' 

-0.03 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

n  =  20,  Q 

-1.48 

-1.06 

-0.55 

-0.19 

0.0 

0.17 

0.20 

Q' 

-0.02 

0.02 

0.0 

0.0 

0.0 

0.0 

0.0 

n  =  oo,  (3 

-1.84 

-1.17 

-0.60 

-0.21 

0.0 

0.12 

0.12 

Q' 

0.03 

0.03 

0.0 

0.0 

0.0 

0.0 

0.0 

Table  2.  Values  of  Q  and  Q'  for  several  values  of  T/ao. 
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n 

02 

03 

-0.1864 

0.7639 

-0.3219 

10 

0.7594 

-0.5221 

20 

0.7438 

-0.6673 

00 

0.6567 

-0.8820 

Table  3.  Polynomial  coefficents  for  Q  —  T  relation 
The  curves  in  Fig.  3  can  be  closely  approximated  by 

The  values  of  O],  oj  and  03,  obtained  by  least  squares  fitting,  are  listed  in  Table  3  for 
several  n  values.  We  have  explored  several  other  values  of  E/ffo  and  u  and  found  that 
the  effect  on  the  value  of  Q  is  negligible.  The  relations  of  Betegon  and  Hancock  [21]  and 
Wang  [25],  correlating  T  with  near-tip  hoop  stress,  also  can  be  rearranged  into  the  form 
of  Equation  3.2. 

To  facilitate  the  use  of  Equation  3.2  in  the  analysis  of  finite  width  geometries  we  have 
provided  normalized  values  of  the  the  stress  intensity  factor  K,  F{a/W),  and  the  T-stress, 
hr(a/W)  and  E{a/W),  for  a  number  of  crack  specimens  in  Table  4.  These  are  taken  from 
Sham  [35,36]  and  Leevers  and  Radon  [37].  The  tabulated  values  allow  us  to  calculate  Q 
in  these  geometries  under  contained  yielding. 

3.2  J  —  T  and  J  —  Q  Approaches 

Two  approaches  to  specifying  families  of  Mode  I  plane  strain  elastic-plastic  crack  tip 
fields  have  been  proposed.  The  first  approach,  suggested  by  Hancock  and  co-workers,  uti¬ 
lizes  the  elastic  T— stress  associated  with  the  second  term  of  Williams’  expansion  (Eq  2.1) 
to  correlate  the  elastic-plastic  crack  tip  fields  of  varying  stress  triaxiality.  The  essence 
of  their  proposal  is  this  —  the  near-tip  stress  triaxiality  Q  can  be  estimated  by  T  using 
Equation  3.2  for  loadings  up  to  fully  yielded  conditions.  We  propose  to  quantify  near  tip 
constraint  using  the  J  —  Q  theory  which  has  a  strong  theoretical  basis  as  discussed  earlier. 

Within  the  MBL  formulation  a  description  of  near-tip  states  by  J  and  Q  is  equi\’a]ent 
to  that  phrzised  in  terms  of  K  smd  T  since  Q  and  T  are  related  by  Equation  3.1  and  J  and 
K  are  related  through 

J  =  (3.3) 
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F(aAV)=K/o°®V^^^ 


CCP 

1.001 

1.006 

1.025 

1.0S8 

1.109 

1.186 

1.303 

1.488 

1.816 

2.581 


DECP 

1.122 

1.121 

1.118 

1.120 

1.132 

1.163 

1.226 

1.343 

1 


TPBB 

1.020 

1.023 

1.027 

1.094 

1.229 

1.465 

1.879 

2.688 

4.627 

12.358 


hT(a/W)  =  T/o“* 


-0.548 

-0.558 

■0.568 


-0.599 


TPBB  t35 


-0.438 

-0.369 

-0.239 

-0.099 

0.075 

0J18 

0.712 

1.487 

3.636 

15.167 


E(a/W)  =  T^/K 


a/W 

CCP  1361 

-1.000  1 
-1.000 

■0.469 

■0.469 

■0.430 

-0.361 

0.20 

-1.004 

■0.473 

-0.233 

0.30 

-1.012 

■0.479 

■0.090 

0.40 

-1.029 

■0.484 

0.061 

0.50 

-1.059 

-0.480 

0217 

0.60 

0.70 

— 

-0.463 

-0.432 

0.379 

0.553 

0.80 

-1.409 

■0.376 

0.786 

0.90 

-1.894 

■0.283 

1227 

Table  4.  Values  of  K,  T  and  S  for  CCP,  DECP  and  TPBB  geometries. 
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under  plane  strain  conditions.  However,  the  equivalence  of  the  two  approaches  does  not 
hold  under  fully  yielded  conditions.  The  J  —  Q  fields  can  exist  over  the  entire  range  of 
plastic  yielding  and  do  not  depend  on  the  existence  of  the  elastic  field  (Eq  2.1).  By  contrast, 
T  is  undefined  under  fully  yielded  conditions.  Furthermore,  our  numerical  investigations, 
to  be  presented  in  the  section  on  finite  width  geometries,  show  that  the  J  —  T  approach 
overestimates  the  actual  stress  triajciality  for  some  geometries  and  imderestimates  it  in 
other  cases  so  that  there  is  not  a  consistent  trend.  Stated  in  another  way,  a  T— stress 
fracture  methodology  could  be  conservative  for  some  geometries  and  non-conservative  in 
others — this  suggests  that  such  an  approach  may  be  impractical. 


Figure  4.  Fracture  specimen  geometries:  (a)  center-crack  panel  (b)  double-edge  crack  panel 
(c)  three-point-bend  bar. 
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4.  EVOLUTION  OF  Q  IN  FINITE  WIDTH  GEOMETRIES 


Q— solutions  have  been  obtained  by  finite  element  calculations  for  the  crack  geometries 
shown  in  Fig.  4.  The  solutions  reported  here  were  obtained  by  anzdyses  based  on  a 
smadl  strain  J2  flow  theory.  A  typical  mesh  for  the  finite- width  crack  geometry  has  about 
1000  four-node  isoparametric  elements.  Readers  are  referred  to  Refs  13  and  14  for  details 
on  finite  element  meshes,  plasticity  formulation  and  solution  procedure.  The  Q  values 
presented  in  the  subsequent  sections  are  evaluated  at  r/(  J/(To)  =  2  unless  stated  otherwise. 
J  is  e\-aluated  using  the  domain  integral  technique  as  discussed  in  Ref  38. 

4.1  Center-Cracked  Panel  (CCP) 

The  evolution  of  Q  for  n  =  3  and  5  and  n  =  10  and  20  is  shown  in  Figs  5  and  6, 
respectively.  Solutions  are  presented  for  the  full  range  of  crack  length  to  width  ratios; 
0.05  <  a/W  <  0.8.  In  the  figures,  J  is  normalized  by  the  crack  length  a  when  a/W  <  0.5 
and  by  the  remaining  ligament  b  when  a/W  >  0.5.  Observe  that  the  stress  triaxiality 
decreases  steadily  with  increasing  J  and  approaches  a  steady-state  slope  at  fully  yielded 
conditions. 

Figure  7a  and  7c  show  the  effect  of  strain  hardening  on  Q  for  a  short  crack,  and  a  deep 
crack,  respectively.  For  both  geometries,  the  loss  of  stress  triaxiality  is  greater  in  the  lower 
hardening  materials. 

The  variation  of  Q  with  distance  is  shown  in  Fig.  7b  eind  7d.  Here,  Q  is  evaluated  at 
r/(J/cro)  =  1,  2,  3 , 4  and  5.  It  can  be  seen  that  Q  has  only  a  slight  dependence  on  r  under 
fully  yielded  conditions.  For  the  range  of  loading  shown  in  Fig.  7,  |Q'|  <  0.03  indicating 
that  the  J  and  Q  are  accurate  descriptors  of  the  field  over  distances  1  <  r/{J/a-Q)  <  5. 

In  Fig.  7b  and  7d  we  also  provide  a  comparison  between  the  zxtual  stress  triaxiality 
and  the  prediction  by  the  T-stress  via  Equation  3.2  The  open  circles  in  Fig.  7b  and 
7d  are  the  T  stress  predictions  and  the  solid  lines  are  the  actual  near-tip  triaxiality.  At 
low  loads.  Equation  3.2  predicts  the  evolution  of  near-tip  stress  triaxiality  accurately. 
However,  at  fully  yielded  conditions,  the  stress  triaxiality  is  incorrectly  predicted.  In  the 
case  of  a/W  =  0.1,  T  underestimates  the  stress  triaxiality  by  about  0.5ao.  For  a  deep 
crack  a/W  =  0.8,  T  overestimates  the  stress  triaxiality  by  a  similar  amount. 
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(a) 


Q 


(c) 


Figure  5.  Center- cracked  pauiel  — evolution  of  Q  with  increasing  J.  a/W  =  0.05,  0.1,  0.2, 


0.3,  0.4  and  0.5,  (a)  n  =  3  (b)  n  =  J  normalized  by  crack  length.  ajW  =  0.5, 


-3.0  -2A  -^0  -15  -1.0  -3.0  -25  -^0  -15  -1.0 

log(J/(a(To))  log(J/(aao)) 


-3.0  -25  -2,0  -15  -1.0  -3.0  -25  -^0  -15  -1.0 

log(J/(bao))  log(J/(b(7o)) 


Figure  6.  Center-cracked  panel  — evolution  of  Q  with  increasing  J.  a/W  =  0.05,  0.1,  0.2, 
0.3,  0.4  and  0.5,  (a)  n  =  10  (b)  n  =  20',  J  normalized  by  crack  length.  a/W  =  0.5, 
0.6,  0.7  and  0.8,  (c)  n  =  3  (d)  n  =  5',  J  normalized  by  remaining  ligament. 


log(J/(aCTo))  log(J/(a(ro)) 


log(J/(baQ))  log(J/(baQ)) 


Figure  7.  Center- cracked  panel.  Effect  of  n  on  the  evolution  of  Q;  (a)  short  crack,  (c)  deep 
crack.  Q  evaluated  at  r^J/co)  =  1,  2,  3,  4  and  5  for  n  =  10;  (b)  short  crack  (d) 
deep  crack.  The  open  circles  are  predictions  baised  on  the  T— stress. 
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4.2  Three-Point  Bend  Bar  (TPBB) 


Solutions  for  Q  for  the  three-point  bend  bar  are  shown  in  Figs.  8  and  9  for  0.05  < 
a/W  <  0.8.  The  behavior  of  Q  in  shallow  cracked  specimens,  a/W  <  0.3,  is  similar  to  that 
seen  for  the  center-cracked-panel,  that  is,  the  loss  of  stress  triaxiality  occurs  gradually. 
When  the  crack  is  sufficiently  deep,  a/W  >  0.3,  high  stress  triaxiality  is  maintained  for 
deformations  characterized  by  J/(aao),  or  J/(bcTo),  less  than  about  0.01.  At  higher  J 
levels,  the  global  bending  stress  field  impinges  on  the  near-tip  region,  r  w  2J/<7o,  causing 
a  rapid  loss  of  stress  triaxiality.  This  occtirs  at  about  J/{b(To)  =  0.02  corresponding  to  a 
deformation  level  which  is  less  than  the  ASTM  limit  for  a  valid  Jjc  test  {J/{b(To)  =  0.04). 

Strain  hardening  effects  on  the  evolution  of  Q  are  displayed  in  Fig.  10a,  10c  and  lOe. 
It  can  be  seen  that  the  effect  of  strain  hardening  on  Q  is  weak  for  deeply  cracked  bend 
bars,  ajW  >  0.4.  The  actual  Q  \^ues  and  the  T— stress  predictions  are  compared  in  Fig. 
10b,  lOd  and  lOf.  It  can  be  seen  that  T  correctly  estimates  the  stress  triaxiality  for  the 
short  crack  geometry  (a/W  =  0.1)  but  fails  to  predict  the  stress  triaxiality  under  large 
scale  yielding  in  the  long  crack  geometries. 

The  behavior  of  Q{f)  over  distances  1  <  r/(J/(To)  <  5  is  shown  in  Fig.  11  for  severed 
deformation  levels  as  measured  by  J/aoQ  or  J/bao.  For  the  short  crack  geometry,  Q  is 
effectively  independent  of  r  over  the  range  of  loading  considered.  In  the  case  of  a/W  =  0.4, 
Q  exhibits  a  dependence  on  r  when  J/aao  >  0.03  (|Q'|  >  0.06);  for  the  deepest  crack 
geometry  this  behavior  occurs  when  J/bcro  >  0.017  (|Q'|  >  0.07).  In  both  cases,  the 
\'ariation  of  Q  with  r  is  linear. 

4.3  Double-Edge  Cracked  Panel  (DECP) 

^-solutions  for  the  double-edge  cracked  panel  are  shown  in  Figs.  12  and  13  for  0.1  < 
a/W  <  0.9.  Surprisingly,  high  stress  triaxiality  imder  fully  yielded  conditions  is  observed 
only  in  the  most  deeply  cracked  geometry,  a/W  =  0.9.  Also  note  that  for  a/W  =  0.9  the 
constraint  falls  slighty  and  then  rises  again  when  the  ligament  is  fully  yielded.  For  shallow 
flaws,  a/W  <  0.5,  Q  is  almost  independent  of  relative  crack  size  a/W. 

Strain  hardening  effects  on  the  evolution  of  Q  are  presented  in  Fig.  14a  and  14c  — 
strain  hardening  effects  are  almost  non-existent  for  a/W  =  0.9.  The  comparison  between 
the  stress  triaxiality  and  the  T-stress  prediction  for  the  double-edged  cracked  pzmel  is 
shown  in  Fig.  14b  and  14d.  The  T— stress  approach  correctly  predicts  the  stress  triaxiality 
for  the  deepest  cracked  geometry  a/W  =  0.9  (see  Fig.  14d)  but  overestimates  the  stress 
triaxiality  for  a/W  =  0.5  (see  Fig.  14b).  Figure  15  shows  the  behavior  of  Q{f)  over 
distzmces  1  <  r/(J/ao)  <  5  for  several  deformation  levels  as  measured  by  J/bco-  It  can  be 
seen  that  Q  exhibits  only  slight  dependence  on  r. 
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Figure  8.  Three-point  bend  bar  — evolution  of  Q  with  increasing  J.  a/W  =  0.05,  0.1  and 
0.2,  (a)  n  =  3  (b)  n  =  5.  a/W  =  0.3  and  0.4,  (c)  n  =  3,  (d)  n  =  5.  a/W  =  0.5, 
0.6,  0.7  £ind  0.8,  (e)  n  =  3  (f)  n  =  5. 
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log(J/(bao))  log(J/(bao)) 

Figure  9.  Three-point  bend  bar  — evolution  of  Q  with  increasing  J.  aJW  =  0.05,  0.1  zmd 
0.2,  (a)  n  =  10  (b)  n  =  20.  afW  =  0.3  and  0.4,  (c)  n  s=  10  (d)  n  =  20.  afW  =  0.5, 
0.6,  0.7  2ind  0.8,  (e)  n  =  10  (f)  n  =  20. 
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Figure  10.  Three-point  bend  bar.  Effect  of  n  on  the  evolution  of  Q;  (a)  short  crack,  (c) 
intermediate  crack,  (e)  deep  crack.  Comparison  of  Q  values  with  predictions 
based  on  the  T— stress;  (b)  short  crack,  (d)  intermediate  crack,  (f)  deep  crack. 
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Figure  11.  Three-point  bend  bar.  Q  evaluated  at  different  positions  ahead  of  crack  tip  for 
several  deformation  levels  as  measured  by  JfLao.  (a)  short  crack,  L  =  a,  (b) 
intermediate  crack,  L  =  a,  (c)  deep  crack  L  =  b. 


log(J/(b<yo))  log(J/(bao)) 


Figure  12.  Double-edge  cracked  peuiel  — evolution  of  Q  with  increasing  J.  a/W  =  0.1,  0.2, 
0.3,  0.4  and  0.5,  (a)  n  =  3  (b)  n  =  5]  J  normedized  by  crack  length.  a/W  =  0.5, 
0.6,  0.7,  0.8,  and  0.9,  (c)  n  =  3  (d)  n  =  5;  J  normalized  by  remaining  ligament. 
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log(J/(ba^)  log(J/(ba^) 


Figure  13.  Double-edge  cracked  panel  — evolution  of  Q  with  increasing  J.  a/W  =  0.1,  0.2, 
0.3,  0.4  and  0.5,  (a)  n  =  10  (b)  n  =  20;  J  normalized  by  crack  length.  a/W  =  0.5, 
0.6,  0.7,  0.8,  and  0.9,  (c)  n  =  10  (d)  n  =  20;  J  normalized  by  remaining  ligament 
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logCJ/CbaJ) 


Figure  14.  Double-edge  cracked  panel.  Effect  of  n  on  the  evolution  of  Q]  (a)  intermediate 
crack,  (c)  deep  crack.  Comparison  of  Q  values  with  predictions  based  on  the 
T-stress;  (b)  intermediate  crack,  (d)  deep  crack. 
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Figure  15.  Double-edge  cracked  panel.  Q  evaluated  at  different  positions  ahead  of  crack  tip 
for  several  deformation  levels  as  meeisured  by  J/6<to.  (a)  intermediate  crack,  (b) 
deep  crack. 


4.4  Use  of  J-Q  solutions 

The  J  —  Q  solutions  provided  in  this  and  the  previous  section  au-e  determined  from 
a  plane  strain  ainadysis.  Therefore  they  are  applicable  to  crack  geometries  which  are  suffi¬ 
ciently  thick  relative  to  the  crack  size  or  uncracked  ligeunent.  In  a  typical  fracture  experi¬ 
ment  J  is  evaluated  from  the  load  displacement  record  and  represents  a  thickness-average 
value.  The  corresponding  value  of  Q  can  then  be  obtained  from  the  figures  presented  or, 
in  the  case  of  well-contained  yielding,  from  Equation  3.2  in  conjimction  with  Table  4.  If 
the  specimen  conditions  at  fracture  do  not  correspond  closely  enough  to  pleine  strain  con¬ 
ditions,  the  Q  solutions  presented  here  are  not  strictly  valid.  The  procedure  may  still  be 
applied  if  we  are  compziring  specimens  of  similiar  thicknesses.  However,  if  test  data  from 
specimens  of  widely  different  thicknesses  are  compared  then  3-D  solutions  for  Q  should  be 
used.  This  aspect  is  discussed  in  Refs  15  and  40  where  a  pointwise  value  of  Q  is  used  to 
quantify  stress  triaxiality  near  a  three- dimensionaJ  crack  front. 

5.  CLEAVAGE  TOUGHNESS  LOCUS 

Kirk  et  al.  [30]  have  measured  cleavage  fracture  toughness  for  A515  steels  at  room 
temperature  over  a  broad  range  of  crack-tip  constraints.  They  tested  edge-cracked  bend 
bars  with  thicknesses  B  =  10,  25.4  and  50.8  mm  and  various  crack  length  to  width  ratios. 
The  measured  toughness  data  is  plotted  against  Q  in  Fig.  16. 

Constraint  effects  on  fracture  toughness  cein  be  predicted  by  using  the  J  —  Q  fields  in 
conjunction  with  a  freicture  criterion  based  on  the  attainment  of  a  critical  stress,  022  = 
at  a  characteristic  microstructural  distance,  r  =  Tc  (Ritchie  c<  al.  [34]).  Within  the  J  —  Q 
annulus  the  normal  stress  ahead  of  the  crack  is  given  by  Equation  2.12  or  more  accurately 
by  Equation  2.13.  In  the  interest  of  simplicity,  we  work  with  the  approximate  form  in 
Equation  2.12.  Suppose  that  Tc  is  within  the  J  —  Q  annulus.  Now  impose  the  RKR 
fracture  criterion  to  get 

-=( - ^  a22{e  =  0)+Q.  (5.1) 

Therefore,  we  can  solve  for  Jc  as  a  function  of  Q  for  selected  \'alues  of  ctc  and  Tc-  Now 
designate  the  toughness  value  at  Q  =  0  as  Jq  and  use  Equation  5.1  to  arrive  at 


Observe  that  the  ratio,  Jc/Jo,  does  not  depend  on  Tc-  The  predicted  variation  of  Jc  with 
Q,  for  (Tc  =  3.5<to,  Jq  =  40  kPa-m  and  n  =  5,  is  indicated  by  the  dashed  line  in  Fig. 
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16.  It  C2in  be  seen  that  the  predicted  toughness  curve  correctly  captures  the  trend  of  the 
experimental  data. 

Sumpter  and  Forbes  [2S]  have  conducted  extensive  testing  on  mild  steel  at  —50®  where 
fracure  occurs  by  cleavage  with  no  prior  stable  crack  growth.  Their  test  program  includes 
shallow  and  deeply  cracked  bend  specimens  and  moderate  to  deeply  cracked  center- cracked 
panels.  Their  results  show  that  the  data  from  fully  yielded  center-cracked  specimens 
are  not  consistent  with  a  J  —  T  toughness  locus  constructed  from  the  bend  specimens; 
this  is  not  unexpected  in  light  of  our  observations  that  the  T-stress  does  not  correctly 
predict  crack- tip  constraint  under  fully  yielded  conditions.  In  contrast,  with  allowance  for 
experimental  scatter,  the  data  from  both  center-cracked  and  bend  specimens  form  a  single 
J  —  Q  toughness  locus. 
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Figure  16.  Cleavage  toughness  data  for  ASTM  A515  Grade  70  steels  tested  at  20® C  using 
edge-cracked  bend  bars;  (Kirk  ei  al.,  [30]).  -I-  for  thickness  B=10  mm,  O  for 
B=25.4  mm,  A  for  B=50.8  mm.  Toughness  curve  predicted  by  J  —  Q  theory  is 
indicated  by  the  dashed  line. 
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6.  J-Q  METHODOLOGY 


The  competition  between  fracture  by  clea\’age  and  ductile  tearing  controls  the  fracture 
resistance  of  ferritic  steels  in  the  ductile-to-brittle  transition  region.  It  is  well  accepted  that 
cleavage  fracture  is  controlled  by  a  critical  hoop  stress  (opening  stress)  at  a  microstruc- 
turally  significant  distance  (Refs  29,  34  and  40).  There  is  cilso  general  agreeement  that  the 
mean  stress  drives  cavity  growth  in  ductile  tearing  (Refs  11  and  24).  Since  Q  quantifies 
both  the  hoop  and  the  mean  stress  relative  to  a  reference  stress  state  it  provides  a  common 
scale  to  interpret  brittle  and  ductile  fracture,  therefore  allowing  both  failure  modes  to  be 
incorporated  in  a  single  toughness  locus. 

Suppose  test  conditions  are  such  that  ductile  and  brittle  mechzmisms  are  operative. 
Fracture  by  cleavage  generally  occurs  at  high  crack  tip  constraint  while  ductile  tearing 
develops  at  low  constraint;  this  is  illustrated  by  the  two  distinct  segments  to  the  toughness 
loci  shown  in  Fig.  17a.  Since  measured  toughness  values  generally  exhibit  scatter,  both  the 
lower  and  higher  toughness  loci  are  indicated  which  define  bounds  for  brittle  and  ductile 
fziilure.  Toughness  values  over  the  full  range  of  crack  tip  constraints  can  be  measured 
by  using  the  test  geometries  depicted  in  Fig.  17a.  As  an  example,  deeply  cracked  bend 
specimens  generate  high  crack  tip  stress  triaxiality,  i.e.  Q  ss  0.  They  produce  driving  force 
curves  which  rise  steeply  and  therefore  intersect  the  toughness  loci  within  a  well-defined, 
narrow  zone  of  the  J  —  Q  diagram.  In  contrast,  center  cracked  panels  md  single  edge 
crack  panels  loaded  in  tension  are  low  constraint  crack  geometries.  They  produce  driving 
force  curves  which  rise  with  shallower  slopes  and  thus  intersect  the  toughness  loci  over  a 
broad  zone  in  the  J  —  Q  diagram.  The  shallow  driving  force  curves  of  low  constraint  crack 
geometry  can  explain  the  considerable  scatter  in  cleavage  values,  Jc,  observed  in  testing 
shallow  crack  specimens  loaded  in  tension. 

Utilization  of  the  toughness  locus  in  fracture  assessments  is  illustrated  in  Fig.  17b. 
Suppose  that  the  material’s  fracture  resistance  under  service  conditions  is  characterized 
by  the  indicated  cleavage- ductile  failure  band.  The  driving  force  curve  for  a  structure  with 
high  crack-tip  constraint,  structure  A,  rises  rapidly  in  the  J  —  Q  space  so  that  cleavage 
fracture  occurs  when  the  driving  force  curve  intersects  the  failure  locus.  In  contrast,  a  low 
constraint  geometry,  structure  B,  induces  a  gradually  rising  driving  force  curve  so  that 
ductile  tearing  is  the  likely  event  at  accidenteJ  overload. 
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Figure  17.  Illustration  oi  J-Q  methodology,  (a)  Laboratory  testing  of  specimens  of  varying 
constraints  to  measure  the  material’s  fracture  resistance.  Circles  indicate  antic¬ 
ipated  scatter  in  the  measured  cleavage  toughness  data  which  define  the  upper 
and  lower  bounds,  (b)  Evaluation  of  structures  using  measured  toughness  locus 
and  predicted  J-Q  load  path  for  two  structural  geometries. 
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